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Abstract 

We study the smoothness of the density of a semiUnear heat equa- 
tion with multiphcative spacetime white noise. Using Malhavin cal- 
culus, we reduce the problem to a question of negative moments of 
solutions of a linear heat equation with multiplicative white noise. 
Then we settle this question by proving that solutions to the linear 
equation have negative moments of all orders. 
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1 Introduction 



Consider the one-dimensional stochastic heat equation on [0, 1] with Dirichlet 
boundary conditions, driven by a two-parameter white noise, and with initial 
condition uq: 

Assume that the coefficients b{t, x, u), a{t, x, u) have linear growth in t, x and 
are Lipschitz functions of m, uniformly in (t, x). 

In [5J Pardoux and Zhang proved that u{t, x) has an absolutely continuous 
distribution for all (t, x) such that t > and x G (0, 1), if cr(0, Uq, Mo(?/o)) 7^ 
for some G (0, 1). Bally and Pardoux have studied the regularity of the 
law of the solution of Equation ([T]) with Neumann boundary conditions on 
[0, 1], assuming that the coefficients h{u) and a{u) are infinitely differentiable 
functions, which are bounded together with their derivatives. 

Let u{t,x) be the solution of Equation ([1]) with Dirichlet boundary con- 
ditions on [0, 1] and assume that the coefficients b and a are infinitely dif- 
ferentiable functions of the variable u with bounded derivatives. The aim of 
this paper is to show that if (t(0, |/o, Mo(|/o)) 7^ for some yo G (0,1), then 
u{t,x) has a smooth density for all {t,x) such that t > and x G (0,1). 
Notice that this is exactly the same nondegeneracy condition imposed in |5] 
to establish the absolute continuity. In order to show this result we make use 
of a general theorem on the existence of negative moments for the solution 
of Equation ([T]) in the case b{t,x,u) = B(t,x)u and a{t,x,u) = H(t,x)u, 
where B and H are some bounded and adapted random fields. 



2 Preliminaries 

First we define white noise W. Let 

W = {W{A), A a Borel subset of \A\ < oo} 

be a Gaussian family of random variables with zero mean and covariance 

E[W{A)W{B)] = \AnB\, 

where |y4| denotes the Lebesgue measure of a Borel subset of M^, defined on 
a complete probability space {Q,J^,P). Then W{t,x) = W{[0,t] x [0,a;]) 
defines a two-parameter Wiener process on [0, cxd)^. 
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We are interested in the following one- dimensional heat equation on [0, oo) x 
[0,1] 

du d'^u d'^W 

with initial condition u{0,x) = uo{x), and Dirichlet boundary conditions 
u(t, 0) = u(t, 1) = 0. We will assume that uq is a continuous function which 
satisfies the boundary conditions uo{0) = mo(1) = 0. This equation is formal 
because the partial derivative |^ does not exist, and ([2]) is usually replaced 
by the evolution equation 

u{t,x) = / Gt{x,y)uo{y)dy+ / / Gt^s{x,y)b{s,y,u{s,y))u{s,y)dyds 
Jo Jo Jo 

+ / Gt-s{x,y)a{s,t,u{s,y))u{s,y)W{dy,ds), (3) 
Jo Jo 

where Gt{x, y) is the fundamental solution of the heat equation on [0, 1] with 
Dirichlet boundary conditions. Equation ([3]) is called the mild form of the 
equation. 

If the coefficients h and a are have linear growth and are Lipschitz func- 
tions of u, uniformly in (t, x), there exists a unique solution of Equation ([3]) 
(see Walsh [8]). 

The Malliavin calculus is an infinite dimensional calculus on a Gaussian 
space, which is mainly applied to establish the regularity of the law of nonlin- 
ear functionals of the underlying Gaussian process. We will briefiy describe 
the basic criteria for existence and smoothness of densities, and we refer to 
Nualart [3j for a more complete presentation of this subject. 

Let S denote the class of smooth random variables of the the form 

F = /(iy(Ai),...,IV(A„)), (4) 

where / belongs to G'^(W^) (/ and all its partial derivatives have polynomial 
growth order), and Ai, . . . , y4„ are Borel subsets of with finite Lebesgue 
measure. The derivative of F is the two-parameter stochastic process defined 
by 

= E • • • ' W^(^n))U(t, X). 

i=i 

In a similar way we define the iterated derivative D^'^^F. The derivative 
operator D (resp. its iteration D^'^^) is a closed operator from into 
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LP{fl;L'^{R'^)) (resp. Lp{Q; L^{R'^^))) for any p > 1. For any p > 1 and for 
any positive integer k we denote by IF''^ the completion of S with respect to 
the norm 



\k,p 



E{\Fn + Y.E 



{Dzi ■ ■■D^.fY dzi ■■■dzj 



Set = DkJ 

Suppose that F = {F\ . . . , F'^) is a d-dimensional random vector whose 
components are in D^'^. Then, we define the MaUiavin matrix of F as the 
random symmetric nonnegative definite matrix 



l<i,j<d 



The basic criteria for the existence and regularity of the density are the 
following: 

Theorem 1 Suppose that F — . . . , F'^) is a d-dimensional random vec- 
tor whose components are in D^'^. Then, 

1. Ifdetap > almost surely, the law of F is absolutely continuous. 

2. IfF' G D°° for each i = 1, . . . ,d and E [(det apy^] < oo for all p>l, 
then the F has an infinitely differentiahle density. 



3 Negative moments 

Theorem 2 Let u{t,x) be the solution to the stochastic heat equation 



dt dx^ ~'~ dtdx ' 
it(0,x) = Uq{x) 

on X & [0,1] with Dirichlet boundary conditions. Assume that B = B{t,x) 
and H — H{t,x) are bounded and adapted processes. Suppose that uo{x) is 
a nonnegative continuous function not identically zero. Then, 

E[u{t,x)-P] < oo 

for all p >2, t > and < x < 1. 



4 



For the proof of this theorem we will make use of the following large 
deviations lemma, which follows from Proposition A. 2, page 530, of Sowers 

Lemma 3 Let w{t,x) be an adapted stochastic process, bounded in absolute 
value by a constant M. Let e > 0. Then, there exist constants Co, Ci > 
such that for all X > and allT > 



P [ sup sup 

vO<t<TO<2:<l 



^0 



Gt~s{x, y)w{s, y)W{ds, dy) 



> A ) < Co exp 



We also need a comparison theorem such as Corollary 2.4 of |6]; see also 
Theorem 3.1 of Mueller [3] or Theorem 2.1 of Donati-Martin and Pardoux 
|2]. Shiga's result is for x G M, but it can easily be extended to the following 
lemma, which deals with x G [0, 1] and Dirichlet boundary conditions. 

Lemma 4 Let Ui{t,x) : i = 1,2 be two solutions of 

dt dx^ * ' ^ dtdx^ 

Ui{0,x) = Uq\x) 

where Bi{t,x), H{t, ) satisfy the same conditions as in Theorem [B 

Also assume that with probability one for all t > 0,x E [0, 1] 



Bi{t,x) < B2{t,x) 

(1)/ X ^ (2)/ N 

Uq [X) < Uq [X). 

Then with probability 1, for all t > 0,x & [0, 1]. 

Ui(t, x) < U2(t, x). 

Proof of Theorem [21 We shall repeatedly use the comparison lemma, 
Lemma m along with the following argument. Observe that if < w{t, x) < 
u{t, x) with probability one, and p > 0, then 

E [u{t,x)-P] < E [w{t,x)-P] . 

Thus, to bound E[u(t, x)~p], it suffices to find a nonnegative function w(t, x) < 
u{t, x) and to prove a bound for E[w{t, x)~^]. Such a function w{t, x) might 
be found using the comparison lemma. Lemma HI 
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Suppose that \B(t,x)\ < K almost surely for some constant K > By 
the comparison lemma, Lemma HI it suffices to consider the solution to the 
equation 

dw d'^w ^ ^ ^ d'^W 

dt dx^ dtdx 

w{0,x) = uo{x) 

on X G [0, 1] with Dirichlet boundary conditions. Indeed, the comparison 
lemma implies that a solution w{t,x) of (171) will be less than or equal to 
a solution u(t,x) of ([5]). Then we can use the argument outlined in the 
previous paragraph to conclude that the boundedness of E [w{t, implies 
the boundedness of E [u{t,x)~P]. 

Set u{t,x) = e~^^w{t,x), where u{t,x) is not the same as earlier in the 
paper. Simple calculus shows that u(t,x) satisfies 

du d'^u ^ d'^W 

dt dx'^ dtdx 
u{0,x) = uo{x) 

and we have 

E [w{t,x)-f] = e^'^E [u{t,x)-P] . 

So, we can assume that K = 0, that is that u{t,x) satisfies ([S]). The mild 
formulation of Equation is 

u{t,x)= Gt{x,y)uo{y)dy + / Gt-s{x,y)H{s,y)u{s,y)W{ds,dy). 
Jo Jo Jo 

Suppose that uo{x) > 5 > for all x G [a, b] C (0, 1). Since ([8]) is linear, 
we may divide this equation by 6, and assume 6=1, and also uo{x) = 
l[a,b](yX). Fix T > 0, and consider a larger interval [a,b] C [c,d] of the form 
d = b + 7T and c = a — 7T, where 7 > 0. We are going to show that 
E{{u{T, x)~^) < 00 for X G [c, d] and for any p > I. Define 

pb+'js 

c = inf inf / Gt{x, y)dy 



0<t+s<T, a--y{t+s)<x<b+-y{t+s) 



a—^s 



and note that < c < 1 for each 7 > and {a,b) G (0,1). Next we 
inductively define a sequence {r„, n > 0} of stopping times and a sequence 
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of processes as follows. Let Vo(t,x) be the solution of (jS]) with initial 

condition uq = l[a,b] and let 

2' 



To = inf < t > : inf Vo{t,x) = - or sup VQ{t,x 

a-7t<z<6+7t 2 0<a;<l C 

Next, assume that we have defined r„_i and for Tn-2 < ^ < Tn-i- 

Then, {f„(t, x), r„_i < t} is defined by ([8]) with initial condition f„(T„„i, x) = 
(f)"l[a-7r„_i,fe+7r„_i](a;). Also, let 



r„ = inf <( t > r„_i : inf Vn{t, x) = ( 

I a—^t<x<b+'yt \ 



or sup Vn(t,x) 

0<x<l 



C\n+l 

2/ 

2 

c 



-71+1 



It is not hard to see that r„ < oo almost surely. Notice that 

/(^N n+l 

inf Vnirn,x) > (-) 
By the comparison lemma, we have that 

u{t, x) > Vn(t, x) 

for all {t, x) and all n > 0. For all p > 1 we have 
E [u{T,xyP] < p(u{T,x) > l) 

n=0 ^ ^ 



(9) 



n+l 



< 



(10) 



Taking into account ([9]), the event {u{T,x) < (f)*^} is included in An 
{r„ < T}. Set <Jn = Tn — Tn-i, for all n > 0, with the convention r_i = 0. 
We have 

sup sup Vi(t, x) > 



P u.< 



n 



Ti-l 



< P 



,r,_i<t<T,_i + ^,0<x<l 



+P ( inf inf Vi{t,x) > (^) 

\Ti^i<t<Ti.i + ^,a-^t<x<b+^t \2J 
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Notice that, for rj_i < t < Tj we have 



Gt-n^^{x,y)dy 



+ 



t rl 



Gt-s{x,y)H{s,y) 



'T-i_i Jo 

As a consequence, by Lemma [3] 



- ) Vi{s,y) 



A - I W{ds,dy). 



P{cr,,<- 
n 



Ti-l 



< P 



sup sup \Ni{t, x)| > 1 

,r,_i<t<r,_i + |,0<x<l 



< Co exp ( —CiW^ 



Next we set up some notation. Let Bn be the event that at least half of 
the variables cXi : i = 0, . . . , n satisfy 



2T 



n 



Note that 



A c E 



since if more than half of the Uj : i = 1, . . . ,n are larger than or equal to 
2T/n then r„ > T. 

For convenience we assume that n = 2k is even, and leave the odd case 
to the reader. Let S„ be all the subsets of {1, . . . , n} of cardinality k = n/2. 
Using Stirling's formula, the reader can verify that as n — ^ oo 



n 
n/2 



0(2") 



:i2) 



Then, 



P{Bn) < P 



< 



i u hf.<v 

\{h,...,ik}es„j=i 
i,...,ifc}eH„ \i=i / 



Using the estimate (fTTi) and (fT^ yields 

P{Bn) < Co2"exp(-Cini/2-.^)- 

< Co exp + Can) 

< Co exp (-Cin=^/2-") 

where the constants Co, Ci may have changed from hne to hne. Hence, 

P (M(r,x) < (0") < Coexp (-Cin3/2--) (13) 
Finally, substituting ([5]) into ([7]) yields E [u{T,x)^'p] < oo. ■ 

4 Smoothness of the density 

Let u{t,x) be the solution to Equation ([2]). Assume that the coefficients b 
and a are continuously differentiable with bounded derivatives. Then u{t,x) 
belongs to the Soboev space D^'^ for all p > I, and the derivative Dg^^uit, x) 
satisfies the following evolution equation 

/■* db 

De,^u{t,x) = J J Gt-s{x,y) — {s,y,u{s,y))De^^u{s,y)dyds 

+a{u{9,0)Gt^9{x,0, (14) 

if 6* < t and Dg^^u{t, x) = if 9 > t. That is, Dg^^u{t, x) is the solution of the 
stochastic partial differential equation 

dDgcu d'^Dg^u db , / nx^^ da, , ^^ -r^ d'^W 
on [^^, oo) X [0,1], with Dirichlet boundary conditions and initial condition 

Theorem 5 Letu{t,x) be the solution of Equation ^ with initial condition 
u{0,x) = uo{x), and Dirichlet boundary conditions u{t,0) = u{t, 1) = 0. We 
will assume that Uq is an a-Holder continuous function for some a > 0, 
which satisfies the boundary conditions Mo(0) = Mo(1) = 0. Assume that 
the coefficients b and a are infinitely differentiable functions with bounded 
derivatives. Then, if a{0,yo,Uo{yo)) ^ for some yo G (0,1), u(t,x) has a 
smooth density for all {t, x) such that t > and x G (0, 1). 
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Proof. From the results proved by Bally and Pardoux in [T] we know 
that u(t,x) belongs to the space for all {t,x). Set 



t,x 



t /.I 



^0 



Then, by Theorem [T] it suffices to show that E{Ci^) < oo for all p> 2. 

Suppose that a{0, yo, Mo(yo)) > 0. By continuity we have that cr(0, y, u{0, y)) > 
(5 > for all y E [a, b] C (0, 1). Then 



f [ {De,^u{t,x)f dide> [ ( [ Dg^^u{t,x)dA d9. 

Ja Jo \Ja J 



Set Yf,^ = Dg^^u{t, x)d^. Fix r < 1 and e > such that e''' < t. From 



Jo 



de + Ct. 



we get 



P{Ct,.<e) < P 



d9 > e 



+P {y,1<V2s'^) 
= P{A^) + P{A2). 

Integrating equation f[T^ in the variable ^ yields the following equation for 
the process {Yf,.,t >9,xe [0, 1]} 

/■* db 

ylx = Gt-six,y) — {s,y,u{s,y))Ylydyds 



+ 



t /•! 



J a 



(15) 



In particular, for 9 = 0, the initial condition is Yq^^ = a{0,^,u{0,^))l[afi]{^), 
and by Theorem [2] the random variable Y^^^ has negative moments of all 
orders. Hence, for all p > 1, 



P{A2) < 



10 



if e < ^o- In order to handle the probabihty P{Ai) we write 



< sup [E 

0<6»<£ 



X 



E 



1/2 



We claim that 



and 



sup E 

o<e<t 



ye 



2g 



< OO, 



sup E 

o<e<e^ 



< e 



2sq 



(16) 



(17) 



for some s > 0. Property (fT6!) follows easily from Equation (fT5|) . On the 
other hand, the difference Yf^ — Y^^ satisfies 



^ t,x ^ t,x 



t 



db 



G,_,(x, y) — is, y, u{s, ?/)) (F^ - Y,"^,)dyds 

t ^1 So- 

Gt.six, y)^{s, y, u{s, y)){Yly - Y^'^,)W{dy, ds) 

Gt-six, y)-g^is, y, uis, y))Y°ydyds 

1 

Gt^six, 2/)^(s, y, u{s, y))Y^ W{dy, ds) 

JO ou 

b 



JO 



+ 
+ 



+ / {a{u{e,0)Gt-e{x,0 - (^{MO)Gt{x,0)d^ 

J a 

5 
i=l 

Applying Gronwall's lemma and standard estimates, to show (1171) it suffices 
to prove that 

sup E{\^,ie)\"')<E''^, (18) 

o<e<e'- 

for z = 3, 4, 5 and for some s > 0. The estimate f|T8l) is clear for z = 3, 4 and 
for i = 5 we use the properties of the heat kernel and the Holder continuity 
of the initial condition uq- Finally, it suffices to choose r > 1 — s and we get 
the desired estimate for P{Ai). The proof is now complete. ■ 
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